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Abstract - The Letter considers dynamics of helical vortices and helical-vortex rings either 
solving directly the equations of motions of a vortex line or using canonical relations following 
from the Hamiltonian equations of motion. An analytical solution in elliptical integrals was found 
for helical-vortex rings in the local-induction approximation. The analysis based on the canonical 
Hamilton relation provides a clear physical explanation for anomalous velocities of helical- vortex 
rings, I.e., for suppression of the velocity and even inversion of its direction at sufficiently large 
amplitude of the helical distortion. The extended local-induction approximation is suggested, 
which provides an exact solution for the equations of motion of helical vortices and helical- vortex 
rings in the limit when the small-pitch helical vortex reduces to a cylindric sheet of uniform 
vorticity. 
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Introduction. — Helical vortices were intensively 
itudied in classical hydrodynamics [Illl]. They appear in 
vakes of propellers and other spinning bodies. A helical 
vortex is a straight vortex line with a circularly polarized 
IjCelvin wave of arbitrary amplitude propagating along it. 
One may roll up the helical vortex into a ring keeping 
ijts helical deformation. Then it will be a helical-vortex 
ring. An ideal vortex ring is an ubiquitous object in vor- 
l^ex dynamics in general and in superfluid turbulence in 
particular and has already been studied about two cen- 
turies. Its dynamics is well established. Recently attention 
was focused on dynamics of helical- vortex rings with large- 
amplitude Kelvin waves propagating around the ring [3H5| . 
This problem being interesting itself may have important 
ipiplications for superfluid turbulence. An interesting out- 
come of studies of this object was that for Kelvin-wave 
amplitudes large enough the vortex ring may move in di- 
rection opposite to its moment. This phenomenon was 
first revealed by Kiknadzc and Mamaladze [3] in numerical 
calculations and the perturbation theory within the local- 
induction approximation. Later it was confirmed \MS\ by 
numerical calculations based on the Bio~Savart law. The 
effect seemed to be mysterious and was called anomalous 
vortex-ring velocity [HE]. 

This Letter addresses dynamics of helical vortices and 
helical-vortex rings in the local-induction approximation 
either solving directly the equations of motions of a vor- 
tex line or using simple canonical relations following from 



the Hamiltonian equations of motion. A fully analytical 
solution in elliptical integrals was found for helical-vortex 
rings with periodic helical distortion of arbitrary ampli- 
tude. It confirms a simplified approach based on the gen- 
eral canonical Hamilton relations. The latter approach 
provides a transparent physical explanation for the origin 
of the anomalous vortex-ring velocity, which now does not 
deserves an adjective "anomalous" . Suppression of the 
ring velocity and the eventual inversion of its direction 
are related with conserved angular momentum around the 
ring axis, which is generated by the helical deformation of 
the vortex ring. Agreement of the canonical approach with 
the full analytical solution (within the area of validity for 
the former) is a proof of its effcctivity. 

The local-induction approximation assumes that a ve- 
locity of a fixed point on a vortex line is induced only by 
the closest segment of the line and is determined by its 
curvature radius. In the case of a helical vortex the clos- 
est segment is restricted by one turn of the helix. The 
presented analysis goes beyond the strict local-induction 
approximation and suggests the extended local-induction 
approximation, which takes into account also the velocity 
induced by other segments (turns) of the vortex line. Ap- 
proximately estimating the effect of the other turns they 
are replaced by a cylindric vortex sheet with uniform dis- 
tribution of vorticity. This approach gives an exact solu- 
tion of the Bio-Savart law of vortex motion in the limit 
of a helix of very low pitch compared to the helix radius. 
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Fig. 1: Helical vortex, (a) Helix of high pitch 2-n /k S> a. (b) 
Helix of low pitch 27r/fc <C a. 



Here Tc is the core radius, is the scale to which the 
velocity field w ~ 1/r penetrates, and = k ln(r,„ /re) /47r 
is the line-tension parameter. The curvature-radius vector 
at any point of the helix is directed along the radius of the 
cylinder while its absolute value is 7^ = (l + fc-^a^)/fc^a. In 
the cylindrical coordinate system (z, r, 0) the unit tangent 
vector has components = 0, = akj \J\ ~V k'^o? , and 
Sz = 1/ a/I + kP-a^ ■ Then solution of the equation ^ of 
motion for the vortex can be presented in two forms (auto 
model solutions): 



or 



(z) = kz — Lot ^ k{z — Vzt) 



(5) 
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Thus the helix motion may be described either as pure 
vertical translation with the velocity Vz, or pure rotation 
with angular velocity ilz- According to eq. ^ 
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Being also accurate in the opposite limit of high pitch, 
where the local-induction approximation is valid, the ex- 
tended local-induction approximation provides a reliable 
interpolation model for a broad spectrum of vortex con- 
figurations. 

Helical vortex. — Let us consider a Kelvin wave of 
arbitrary amplitude propagating along a straight vortex 
line. In the Cartesian coordinate frame with the z axis co- 
inciding with the unperturbed vortex line, the transverse 
displacements in the circularly polarized Kelvin wave are 



a cos{kz — Lot) , y ^ a shi{kz — Lot) . 



(1) 



So in the Kelvin mode the vortex line forms a helix moving 
along a cylinder coaxial with the z axis with the pitch 
equal to the Kelvin wavelength 27r/fc (fig.[l|). 

In the local induction approximation the equation of 
motion for a curved vortex line is 



PsKs{Rl) X 



~dt 



n 



(2) 



Here ps is the supcrfiuid mass density, k = h/m is the 
circulation quantum, the position vector for points on 
the vortex line, s{Rl) is the unit vector tangent to the 
vortex line at the point with the position vector Rl, TZ 
is the curvature-radius vector, which is determine by the 
relation 

^ - ^ (3) 

dl is the element of the vortex line length, and e is the 
vortex line tension equal to the vortex-line energy per unit 
length: 
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In the limit of small a this yields the dispersion relation 
oj = Vgk^ for the linear Kelvin wave. But the frequency 
decreases with increasing amplitude a of the Kelvin wave. 

Expressions ([7]) and 1^ can be also derived from the 
general canonic Hamilton equations 
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where Pz and Mz arc the linear momentum and the an- 
gular momentum along the z axis. 
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is the energy of the vortex line, and L and Rc are the 
hight and the radius of the container with the helical vor- 
tex. Partial derivatives in eq. (O took into account that 
the only varying parameter of the helix is its radius a. 
Indeed, imposing the periodic boundary conditions in the 
box of the hight L the pitch 27r/fc of the helix becomes a 
topological invariant. The first term in the expression (jlOp 
is the kinetic energy at the distances r < y/l + k^a^/k, 
where the helical deformation increases the length of the 
vortex line by the factor + fc^a^, while the second term 
is the kinetic energy at distances r > \/l + k'^a^ jk, where 
helical deformation does not affect the velocity v — nj^-nr 
determined only by the circulation quantum. 

In general a momentum of a vortex line is given by [T] 



^ / [flL X dl\ 



(11) 
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Fig. 2: Helical-vortex ring 



Here dl — sdl, and dl is the element of the vortex line 
length. The expression leaves the vortex momentum un- 
defined (similar to an undefined dipole moment of a sin- 
gle electron). But physically only the difference of the 
momenta of a distorted and a straight vortex along the 
chosen vertical axis is of importance. Then 



(12) 

This momentum is nothing else as the momentum of the 
constant vertical velocity field /tfc/27r inside the cylinder 
of radius a around which the vortex line is winding. The 
winding vortex line with tilted circulation vector produces 
a jump of the vertical velocity equal to hik/2'n in average, 
and the latter is equal to the vertical velocity inside the 
helix under natural assumption that there is no velocity 
outside the helix. Though the vertical velocity determines 
the whole momentum, the kinetic energy of the vertical 
flow is ignored in the local induction approximation since 
it has not a large logarithmic factor. But we shall take 
the kinetic energy of the vertical flow into account study- 
ing the case when the strict local-induction approximation 
becomes invalid (see below). Similarly, calculating the z 
component of the angular momentum of the velocity 
field induced by the helical vortex only the difference be- 
tween the angular momenta of a distorted and a straight 
vortex is of importance, and 



PskL 2 
Mz = — a 



(13) 



Note a negative sign in this expression. Indeed, the an- 
gular momentum of a vortex line displaced from the con- 
tainer axis is less than the angular momentum of a vortex 
line exactly along the container axis. The latter is very 
large and depends on the container radius, but does not 
affect the outcome of the analysis determined only by mo- 
mentum difference. 



Helical- vortex ring: canonical formalism. Now 

let us roll up the helical vortex into a ring of the radius 
R = L/2t:. The helix makes n turns around the vortex 
ring (fig. [2]). It is assumed that the radius a of the he- 
lix is much smaller than R but may be of the same order 
or larger than the pitch 27ri?/n of the helix (the Kelvin 
wavelength). We shall derive the velocity of the ring in 
the absence of dissipation using the canonical relation for 
the ring velocity, which takes into account that the coor- 
dinate along the ring axis and the linear momentum along 
the same axis constitute a pair of canonically conjugate 
Hamiltonian variables: 



dE _ DE/dR 
aP ^ dP/dR 



(14) 



The momentum P = ttPskR^ related with the ring is not 
affected by helical deformation of the vortex line. The 
expression (jlOp for the energy given in the previous section 
can be used now after replacing L by 2t:R and k by n/R. 
So the energy is 
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But now the system has two degrees of freedom, and one 
should decide what is going on with the vortex helix if the 
ring radius varies. Certainly the number n of turns cannot 
vary being a topological invariant. Another invariant is 
the angular momentum along the ring axis. For large ring 
radius i? ^ a it can be estimated from the component P^ 
of the straight helical vortex, eq. (|12p : 
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Variation of the energy with respect to P must be per- 
formed at fixed Max- So the growth of R is accompanied 
by a decrease of the helix radius a: da/dR = —a/2R. 
Keeping the most important logarithmic terms the ring 
velocity is 



K 1 - n^a^/2R^ , ^/WT^n?a? 

Vr ^ , — In 

47r Vi?^ + n^a^ 



(17) 



So at large number of helix turns n > R/a the sign of the 
velocity changes, and the ring starts to move in the direc- 
tion opposite to its momentum. This is a phenomenon of 
the anomalous vortex-ring velocity, which is here directly 
explained by the canonical equations of motion and the an- 
gular momentum conservation law. Expansion of eq. (|17p 
in na/R yields the perturbation-theory result by Kiknadze 
and Mamaladzc The critical value of na/R = \p2 at 
which the vortex ring stops is larger than najR— 1 given 
by the perturbation theory and is in reasonable agreement 
with the numerical calculations with the Bio-Savart law 
giving the critical value na/R = 1.7 at large n [3]. 

The angular velocity of the helical-vortex ring also fol- 
lows from the canonical relation: 
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(18) 
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Now taking the partial derivative the hnear momentum 
P {i.e., the average radius R) must be fixed. Rotation 
with angular velocity 57 corresponds to translation of the 
straight helical vortex, so 57 = Vz/R, where Vz is given by 
eq. ([7]) but with k replaced with n/R. 

Helical-vortex ring: general analytical solution. 

— The problem of a helical-vortex ring in the local in- 
duction approximation can be solved analytically with- 
out any assumption on the ratio a/R. The shape of the 
helical-vortex ring is determined in the cylindrical co- 
ordinates by the two functions r{(p,t) and z{(t),t). We 
look for an automodel solution r{4>,t) — r{(f> — iU) and 
z{(j),t) — z{(j) — nt)+Vrt, which describes stationary trans- 
lation and rotation of the ring with the linear velocity Vr 
and the angular velocity 57 respectively. The vector equa- 
tion ([2]) of motion reduces to the two equations for the two 
time-independent functions r(^) and z((j)): 
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Or- 



dr 



r^l -t- (dr/rd0)2 + (dz/rd(^)2 



0, (19) 



dr 



[1 + (dr/rd(^)2 -h (dz/rd0)2]i/2 

These are two equations of balance for forces along the axis 
z and for moments around the axis z (azimuthal forces) 
applied to an element of the vortex line: The left-hand 
sides present the Magnus forces, whereas the right-hand 
sides are the line-tension forces proportional to the com- 
ponents of the curvature vector in the cylindric coordinate 
frame. The first integration of the equations is straight- 
forward: 

2 



Not all necessary boundary conditions were still satis- 
fied. In order to obtain a continuous curve spiraling along 
toroidal surface one should require that at the maximum 
r ~ i?2 the coordinate z returns to the same value as at 
the minimum r = Ri. The latter chosen to be zero this 
yields the condition z{R2) = 0. The second condition is 
dictated by a chosen period of the helix around the ring 
and is given by (/)(i?2) = 7r/n, where n is the number of 
turns by the helical vortex along the ring. All together 
. (20)this provides a full analytic solution of our problem. 

First one can check the limit of a weak Kelvin wave 
propagating around the ring (the linear theory). In this 
limit a = R — Rl = R — R2 is very small. Then our 
solution yields that the translational velocity Vr = i^s/Roi 
the ring does not differ from that of an unperturbed ring. 
At the same time the angular velocity of the helix is given 

by n 
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d(^; [Vr{r^ - Rl) + i?i^4j.2_f^2(^2_^2)2]2 



(21) 
(i) 



At integration two boundary conditions were used: 
dz/d(j) = at r = i? [the maximum of the function z{<j))\, 
and (ii) dr/dcj) = at r = i?i < R [the minimum of the 
function r{(j))]. The maximum of the function r((/)) (the 
the vortex line position most distant from the ring axis) 
is at the distance r = R2 given by 

5*66 eauation \2'2\ next page 

The second integration yields expressions for functions 
(j){r) and z{r) via elliptic integrals: 

, , , _ Vr{Rl+Rl) 



nR^^/R.: 



[F{K,m) 



n 



(Rl 



+ Rl)U 



Rl' 



VAvi-n-^{Rl- R-^y 
nRiRiVRl + Rl 

- RlF{K,m) 



(23) 



Vs\/v?- — 1/i?^, like that predicted by the linear 
theory of the Kelvin wave along the circular vortex loop. 
The most interesting case for us is when the amplitude 
of the Kelvin wave is still small compared with the ring 
radius but not small compared with the Kelvin wavelength 
27ri?/n. This requires large winding numbers n ^ 1. In 
this limit one obtains exactly the same result as given by 
cq. (jl7p derived from the canonical formalism. In order 
to demonstrate this one should keep terms of the second 
order m. a ~ R — R\ and therefore expand the elliptic 
integrals in cq. (j25p up to terms . For not large wave 
numbers n and large amplitudes a ~ R ~ R\ one should 
use values of unexpanded elliptic integrals in calculations. 
Figure [3] shows helical- vortex rings of wave number n = 4 
for various ratios R\ / R. The vortex line is winding around 
a surface of a toroid, which changes from an ideal torus at 
Rl ^ R with circular cross-section to a toroid obtained 
by revolution of a figure of eight at i?i — > (fig. HJ. 

Our analysis demonstrates the existence of an auto- 
model solution when the whole motion of the helical- 
vortex ring reduces to its pure translation and rotation 
without any variation of its shape. Meanwhile, Barenghi 
et al. [4] reported that in their numerical calculations us- 
ing the Bio-Savart law they detected oscillations of the 
total length of the vortex line with time. It is difficult 
to imagine that the automodel solution exists only in the 
local-induction approximation, especially keeping in mind 
that calculations on the basis of the Bio-Savart law con- 
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Fig. 3: Helical-vortex rings of four-fold symmetry n = 4 for 
various ratios Ri/R. The left column shows projections of 
rings on the plane normal to the ring axis (the axis z). Dashed 
lines show unperturbed ideally circular rings of the radius R. 
The right column shows 3D images of rings. 
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Fig. 4: Axial cross-section of toroids over which vortex lines 
wind for the case of four-fold symmetry n = 4 and ratios 
Ri/R = 0.9 ,0.6, and (shown near the curves). The cross- 
section varies from two ideal circles at 7?i — ^ 7? to the figure 
eight at J?i ^0. The latter corresponds to the bifurcation 
point where reconnection to a few separate loops may occur. 



firm its existence for straight helical vortices Appar- 
ently, in simulations by Barenghi ct al. the vortex line 
length oscillates because of excitation of more than one 
Kelvin mode. 

Beyond the local induction approximation. 

All the results received above were obtained strictly 
within the local-induction approximation, which takes into 
account only the vortex velocity induced by the close seg- 
ment of the vortex line and determined by its curvature.. 
Meanwhile, in a helical vortex of low pitch 27r/fc <C a 
[fig. [Ijb)] distant parts of the same vortex line (distant in 
the sense of the distance measured along the vortex line) 
approach to a given point of the vortex line to distances 
much shorter than the curvature radius. Then the local- 
induction approximation in its strict meaning becomes in- 
valid: one must take into account velocity induced by the 
other turns of the spiral. The most rigorous way to deal 
with such cases is using the Bio-Savart law. In classical 
hydrodynamics an analytical solution of the Bio-Savart 
equations for helical vortices was obtained in the form of 
infinite series of Bessel functions (Kapteyn series) [2] . On 
the other hand, one may suggest a much simpler but still 
reasonably accurate approach, which takes into account 
the velocity induced by distant parts of the vortex line by 
modification or extension of the local-induction approxi- 
mation, 

We start from the straight vortex shown in fig. [Ub). 
At low pitch the surface of the cylinder around which the 
vortex line is winding can be considered as a continuous 
vortex sheet with the jump of the vertical velocity com- 
ponent Kk/2'K. Assuming that there is no velocity outside 
the helix, but only inside, one obtains the kinetic energy 
of axial (vertical) motion per unit length: psK^fc^a^/87r. 
Adding this energy to the energy in the local-induction 
approximation given by eq. (jlOp one obtains 

47r V VYT^ 

(26) 

The added kinetic energy (the last term in the parenthe- 
ses) is related to the velocity induced by all turns of the 
vortex line around the cylinder of radius a. It is beyond 
the strict local-induction approximation since it does not 
contain a large-logarithm factor. But it provides the most 
essential contribution at afc 3> | ln(fcrc)|. In this limit this 
approach must yield the same results as the Bio-Savart 
law. So the expression ([^5]) is a reliable interpolation be- 
tween two limits ak <C |ln(fcrc)| (local-induction limit) 
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and ak ^ |ln(fcrc)| (continuous- vortex-sheet hmit). The 
extended local-induction approximation does not change 
expressions for the linear momentum [eq. ()12[)] and the 
angular momentum [eq. ()13p ] since at their derivation 
the local-induction approximation was not used. 

Using the expression p6[) for the energy in the canonic 
expressions for the translational and angular velocities 
[eq. dH])] one obtains: 



u.k 



Vz = -r 



nk 



An 



(27) 



(28) 



VT+k^ 

Note that the continuous- vortex-sheet contribution nk / An 
to the translational velocity is an average between the zero 
velocity outside and the velocity Kk/2n inside the cylindric 
vortex sheet like in the case of plane vortex sheets well 
known in hydrodynamics. 

A similar modification of the local- induction approxima- 
tion is possible for helical- vortex rings. The added contin- 
uous velocity is now an azimuthal velocity inside the helix 
formed along the ring. This leads to modification of the 
expression eq. (|17p for the ring translational velocity: 



Vr = =- m 

An V-R^ + n^a^ nVc 



87ri?3 



(29) 



A negative sign of the new non-logarithmic term demon- 
strates that the reverse motion of the vortex ring due to its 
helical deformation remains also in the continuous- vortex- 
sheet limit. 

The principle, at which the extended local-induction ap- 
proximation is based, is not new. The well-known Hall- 
Vinen-Bekarcvich-Khalatnikov theory for rotating super- 
fluids [71 IH] was based on the same principle: The log- 
arithmically large self-induction line-tension contribution 
to the vortex-line velocity is combined with the velocity 
induced by other vortices in the vortex array, which is ap- 
proximated by continuous vorticity. In our case there is no 
other vortices, but the continuous vorticity is introduced 
as an approximation for the effect of other turns of the 
same vortex line. 

Other extensions and implications of the analy- 
sis. — Studying the helical- vortex rings the analysis ad- 
dressed the case of a single Kelvin mode with fixed wave 
number. A natural extension of the analysis is the case 
of the ensemble of Kelvin modes assuming that the effect 
of various modes is additive [3], i.e., replacing n^a^ by 
J2 n^al- Then C = 2n fWTY/n^ is the total length 



of the vortex line increased by the vortex line fluctuations. 

We have obtained that elongation of the vortex line de- 
creases the translational ring velocity, but this conclusion 
is not universal being related with conservation of the an- 
gular momentum at propagation of an isolated single vor- 
tex ring. Dynamics of vortex rings is widely exploited in 



the theory of superfluid turbulence. One may not expect 
that in a vortex tangle with many interacting rings the 
conservation law for a particular vortex is applicable. On 
the other hand, a possibility can be considered that inter- 
action of Kelvin modes on various rings establishes some 
fixed average energy per unit length of any vortex line. 
Then the partial derivative of the energy with respect to 
the ring radius must be calculated at fixed ratio ^ n^a^f R 

n 

and the translational velocity Vr = VsC/2nR^ exceeds that 
of the ideal circular vortex ring. One may consider this 
as renormalization of the line tension parameter Vs by the 
factor C/2nR taking into account elongation of the vortex 
line by short-wavelength Kelvin modes. 

Conclusions. — The Letter studies dynamics of heli- 
cal vortices and helical- vortex rings in the local-induction 
approximation either analytically solving directly the 
equations of motions of a vortex line or using canonical 
relations following from the Hamiltonian equations of mo- 
tion. The analytical solution of equations of motions fully 
confirms the results following from simple canonical rela- 
tions. The analysis based on these relation explains the 
origin of the anomalous vortex-ring velocity, which docs 
not deserves an adjective "anomalous" anymore. Sup- 
pression of the ring velocity and eventual inversion of its 
direction are related with conserved angular momentum 
around the ring axis, which is generated by the helical de- 
formation of the vortex ring. The analysis goes beyond the 
local-induction approximation and suggests the extended 
local-induction approximation, which agrees with the ex- 
act solution for helical vortices and helical- vortex rings in 
the limit when the vortex line winding around a cylinder 
(helical vortex) or a torus (helical-vortex ring) reduces to 
a sheet with uniformly distributed vorticity. 
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